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ABSTRACT 

We present an analysis of the first high-resolution spectra measured from an accretion-driven millisecond X- 
ray pulsar in outburst. We observed XTE J1751-305 with XMM-Newton on 2002 April 7 for approximately 
35 ks. Using a simple absorbed blackbody plus power-law model, we measure an unabsorbed flux of 6.6 ± 
0.1 x 10" 10 erg cm" 2 s" 1 (0.5-10.0 keV). A hard power-law component (r = 1.44 ±0.01) contributes 83% of 
the unabsorbed flux in the 0.5-10.0 keV band, but a blackbody component (kT = 1.05 ±0.01 keV) is required. 
We find no clear evidence for narrow or broad emission or absorption lines in the time-averaged spectra, and the 
sensitivity of this observation has allowed us to set constraining upper-limits on the strength of important features. 
The lack of line features is at odds with spectra measured from some other X-ray binaries which share some 
similarities with XTE J1751-305. We discuss the implications of these findings on the accretion flow geometry 
inXTEJ1751-305. 



very low (~ 0.01 M Q ). 

After the discovery of XTE J175 1-305, we submitted a 
target-of-opportunity request to XMM-Newton for the purpose 
of studying this source with high-resolution spectroscopy. The 
X-ray spectrum of the first system — SAX J1808.4-3658 — 
could only be studied during outburst using RXTE and with the 
Wide Field Cameras abomdBeppoSAX, which have only mod- 
erate spectral resolution. In this Letter, we present an analy- 
sis of the time-averaged EPIC-pn and Reflection Grating Spec- 
trometer (RGS) data. The spectra obtained represent the first 
CCD- and gratings-resolution measurements from a millisec- 
ond X-ray pulsar in outburst. 

2. OBSERVATION AND DATA REDUCTION 

XMM-Newton observed XTE J1751-305 beginning on 2002 
April 07.5 (UT). The EPIC-pn camera was operated in "tim- 
ing" mode to prevent photon pile-up effects and to allow pulse- 
phase-resolved spectroscopy . The RGS was operated in the 
standard "spectroscopy" mode. The "thin" optical blocking fil- 
ter was selected for the EPIC cameras. The EPIC-MOS cam- 
eras were not operated in modes optimized for a source of this 
intensity, and we do not consider the MOS spectra here. 

The EPIC-pn and RGS data were reduced using SAS version 
5.3. We applied the standard reduction procedure "epproc" to 
produce a calibrated pn event list. In "timing" mode, the spatial 
information is compressed into one dimension. We extracted 
source data in a rectangle (with a width of 37") along the length 
of the CCD, and background data from adjacent regions; only 
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1. INTRODUCTION 

Millisecond radio pulsars are thought to be created in neu- 
tron star low-mass X-ray binaries (LMXBs). In those LMXBs, 
accreting matter may spin-up the neutron star (see, e.g., Bhat- 
tacharya & van den Heuvel 1991). Therefore, it is expected 
that millisecond pulsars should also be found in LMXBs dur- 
ing the accretion phase of the binary as X-ray pulsars. Although 
evidence for rapidly spinning neutron stars in LMXBs was in- 
ferred from the burst oscillations which were seen during type- 
I X-ray bursts in several systems (see Strohmayer 2001 for a 
review), the detection of millisecond pulsations in persistent 
emission remained elusive for many years. In 1998, the first 
such system was discovered (SAX J1808.4-3658, which has a 
spin frequency of 401 Hz; Wijnands & van der Klis 1998a). 
This source was extensively studied due to its obvious impor- 
tance for binary evolution scenarios and accretion flow geome- 
tries and dynamics (see, e.g., Wijnands & van der Klis 1998b 
and references therein). 

In the spring of 2002, Markwardt & Swank (2002a) reported 
the discovery of the second accretion-driven millisecond pulsar 
XTE J1751-305 (435 Hz) in the Rossi X-ray Timing Explorer 
(RXTE) bulge scan observation program. About one month af- 
ter the discovery of this system, XTE J0929-314 was discov- 
ered (185 Hz; Remillard, Swank, & Strohmayer 2002.). The 
neutron stars in both systems are in orbit around a compan- 
ion star, with an orbital period of ^42 minutes (Markwardt & 
Swank 2002b; Galloway et al. 2002). These systems are very 
tight binaries and the inferred mass of their companion stars is 
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"single" and "double" events were selected. The net pn ex- 
posure was 33.7 ks. Using LHEASOFT version 5.1, the tool 
"grppha" was used to rebin the spectrum to require a minimum 
of 20 counts per bin. We used the trial "timing" mode redistri- 
bution and ancillary response functions developed by the EPIC- 
pn calibration team. We applied the standard reduction proce- 
dure "rgsproc" to produce calibrated RGS event lists, first-order 
spectra, and response functions. Periods of high instrumen- 
tal background were excluded, giving a net RGS exposure of 
34.2 ks. The RGS spectra were grouped to require at least 20 
counts per bin. 

3. ANALYSIS AND RESULTS 

Preliminary spectroscopic results from this XMM-Newton 
observation were reported in Miller et al. 2002a and Miller et al. 
2002b. The accurate position of XTE J1751-305 was reported 
by Ehle et al. (2002), also based on analysis of this observation 
(R.A. = 17/i51ml3.5s,Decl. = -30°37'22", equinox 2000.0). 

The internal calibration accuracy of the EPIC cameras is bet- 
ter than 5% for on-axis sources (Kirsch 2002). We find devia- 
tions up to 0.6 keV which are not removed by spectral model- 
ing, and therefore restrict our analysis of the pn spectrum to the 
0.6-10.0 keV band. The average count rate (79.36 counts s" 1 ) is 
well below that at which photon pile-up occurs in timing mode 
(1500 counts s -1 ). Tests of the effective area calibration of the 
RGS reveal accuracies of 5% in the 7-36A o (0.35-1.77 keV) 
range and 10% at wavelengths shorter than 7 A (energies above 
1.77 keV; den Herder 2002). Although our best-fit model for 
the EPIC-pn spectrum (see below) provides an adequate fit to 
the RGS spectrum (x 2 /^ = 1.314, = 1463), broad deviations 
are apparent at energies above 1.8 keV. We therefore rely on 
the EPIC-pn spectrum to characterize the continuum emission, 
and search the RGS band in contiguous 3A slices for narrow 
features. 

The EPIC-pn and RGS spectra were analyzed using XSPEC 
version 11.1. All errors quoted in this paper are 90% confidence 
errors. Note that systematic errors are not added to the spectra 
to account for flux calibration uncertainties. 

3.1. The EPIC-pn Spectrum 

We first considered a model consisting of blackbody and 
power-law components, modified by neutral absorption in the 
ISM (using the "phabs" model in XSPEC). We measure a 
best-fit equivalent hydrogen column density of Nh = (9.8 ± 
0.1) x 10 21 atoms cm" 2 . The blackbody temperature ob- 
tained is moderate: kT = 1.05 ±0.01 keV. For spherical sym- 
metry, this translates into a blackbody radius of R = f 2 x 
3.01^q^ (d/10 kpc) km, where the spectral hardening fac- 
tor / is the ratio of the color temperature and the effective 
temperature, and d is the distance to the source. The mea- 
sured power-law index is T = 1.44 ±0.01; the normalization 
is (6.4 ±0.1) x 10" 2 (photons keV -1 cm" 2 s" 1 at 1 keV). This 
model produces an adequate fit to the data: x 2 / u = 1.107 
(where v is the number of degrees of freedom; v = 1884 for this 
model; see Figure 1). We note that a simple power-law model 
returns a very poor fit (x 2 /^ > 4, v = 1886); the soft blackbody 
component is strongly required by the data. 

With this simple but standard model, we obtain an unab- 
sorbed 0.5-10.0 keV flux of (6.6 ±0.1) x 10~ 10 erg cm 2 s" 1 , 
or (0. 170 ± 0.004) photons cm" 2 s" 1 . The hard power-law com- 
ponent contributes 83% of the unabsorbed flux 
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FIG. 1.— The EPIC-pn spectrum of XTE J175 1-305, fit with a 
model consisting of blackbody and power-law components modified 
by photoelectric absorption (see Section 3 for details). The spectrum 
and data/model ratio shown above are rebinned for visual clarity. 

in the 0.5-10.0 keV band (an extrapolation from the 0.6- 
10.0 keV band as the 0.5-10.0 keV range is more conventional). 

Motivated by evidence for a Comptonizing volume in 
SAX J1808.4-3658 (Gierlinski, Done, & Barret 2002), we also 
fit the EPIC-pn spectrum with the "compTT" Comptonization 
model (Titarchuk 1994), which attempts to account for the 
Compton-upscattering of cool photons in a hot corona self- 
consistently. This model does not provide an acceptable fit by 
itself (x 2 / v = 2-4, v = 1884). We find a degeneracy between a 
cool corona (kT ~ 3 keV) with high optical depth (r <~ 8), and 
a hot corona (kT ~ 40 keV) with lower optical depth (r ~ 1), 
each up-scattering a seed photon distribution peaking between 
kT =0.5-0.6 keV. Finally, we tried fitting a blackbody compo- 
nent in addition to compTT. This model might correspond to a 
scenario in which the blackbody emission region is only par- 
tially covered by the corona. This model also yielded a poor fit 
(X 2 /V~2). 

The strength, ionization state, and profile of Fe Ka emis- 
sion lines — if produced through irradiation of the accretion 
disk — can reveal a great deal about the accretion geometry of 
low-mass X-ray binaries. Asai et al. (2000) report the possible 
detection of an Fe Ka line in 4U 1820-30 (P or b = 1 1 .4 minutes; 
Stella, White, & Priedhorsky 1987). The line was found to be 
ionized (E = 6.6 ± 0. 1 keV), broad (FWHM = 0.7!j> 2 keV), and 
weak (Wxa = 33±} 2 eV). Assuming a line of equivalent FWHM 
in the 6.40-6.97 keV range (Fe I — Fe XXVI), we measure a 
95% confidence upper-limit of just 6 eV in the EPIC-pn spec- 
trum of XTE J175 1—305. Using a Gaussian with zero width 
to model a line narrower than the EPIC-pn resolution, the 95% 
confidence upper limit on the equivalent width of any such nar- 
row Fe Ka lines is less than 4 eV (6.40-6.97 keV). 

3.2. Pulse-Phase-Resolved Spectroscopy: First Results 

We applied the SAS task "barycen", and using the binary 
system parameters reported by Markwardt et al. (2002) we pro- 
duced spectra from the "high" and "low" parts of the pulse. We 
find that the model applied to the time-averaged spectrum is an 
acceptable description of these spectra (no convincing narrow 
or broad emission or absorption features are found). Interest- 
ingly, the slope of the power-law component is consistent in the 
two spectra but the normalization changes, and the blackbody 
component changes in temperature but not in normalization. 
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FIG. 2.— The RGS spectra of XTE J1751-305 (RGS-1 in black, 
RGS-2 in red; both are rebinned for clarity and narrow features are 
instrumental) fit locally in the region of the Ne K edge (14.25A) with 
a simple power-law plus edge model. The edge depth was fixed at the 
solar value but the power-law index and normalization were allowed to 
vary; a trend in the data/model ratio is clearly visible. Fits made with 
a variable edge depth suggest Ne is likely under- abundant (see text). 

Similar behavior was discovered in pulse-phase-resolved 
spectra of SAX J1808.4-3658 (Gierlinski, Done, & Barret 
2002). 

We note that a possible feature at approximately 2.9 keV ap- 
pears to vary in width, intensity, and centroid energy in the 
"high" and "low" spectra; however, the error limits on the line 
parameters overlap. This feature is also visible in the time- 
averaged spectrum (see Figure 1). As there are no astrophys- 
ically abundant elements with transitions near this energy, it is 
tempting to associate this feature with a red-shifted absorption 
line from the neutron star surface (given the temperature of the 
blackbody, perhaps Fe XXV or Fe XXVI). Although an F-test 
finds that the addition of a Gaussian to model this feature is 
significant at the la level of confidence in the "low" part of 
the pulse (3(7 in the "high" part), a similar feature is apparent 
when the EPIC-pn spectrum of the Galactic black hole candi- 
date XTE J 1650-500 (Miller et al. 2002c) is fit using the same 
response matrix. It is likely that this feature is due to a defect 
in the timing-mode response matrix. 

3.3. The RGS Spectra 

We analyzed the RGS spectra in contiguous 3A slices for 
narrow emission or absorption lines. In each 3 A slice, a simple 
power-law plus ISM edge(s) model was fit to the data. We find 
no convincing evidence for narrow absorption features in this 
band (0.3-2.5 keV, or 5^0 A). At energies below -0.6 keV 
(wavelengths longer than ~21A) the sensitivity is low due to 
the relatively high column density along the line of sight to XTE 
J1751-305. 

Recently, evidence for enhanced Ne abundances has been 
found in compact binaries (Schulz et al. 2001; Juett, Psaltis, & 
Chakrabarty 2001). We fit the RGS spectra around the position 
of the Ne K edge (14. 25 A ) with a power-law model between 
13.25A and 15.25A(see Figure 2). The power-law index and 
normalization were both allowed to vary in this fit. Using the 
cross-sections of Verner et al. (1993) and solar abundances rel- 
ative to H as per Morrison & McCammon (1983), we find that 
Ne may be imc/er-abundant in XTE J1751-305: 
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FIG. 3.— The RGS-2 spectrum of XTE J1751-305 in the region of 
the Ne X Ly-a line at 12.13A, fit with a simple power-law plus Gaus- 
sian model. Our best fit model does not require an emission line within 
12.13 ± 0.1 A, assuming a FWHM equivalent to the prominent Ne X 
Ly-a emission line seen in the 42-minute X-ray pulsar 4U 1626-67 
(Schulz et al. 2001). Narrow features are instrumental. 

the 95% confidence upper-limit is that the abundance of Ne 
is only 77% of the solar value. A solar Ne abundance does not 
allow for an optimal characterization of the local spectrum (see 
Figure 2). 

The binary X-ray pulsar 4U 1 626-67 also has a 42-minute or- 
bital period. Moreover, the distance to this system may be close 
to 8 kpc (Chakrabarty 1998); we anticipate a Galactic center 
location for XTE J1751-305 and therefore a similar distance. 
A C 'hand ra/HETGS spectrum of 4U 1626-67 reveals broad- 
ened emission lines and Doppler-shifted line pairs (Schulz et al. 
2001). The Ne Ly-a line (12.13A) is particularly strong (a flux 
of 8. 15 ±0.93 x 10" 5 photons cm -2 s _1 is measured; FWHM 
= 2860 ± 330 km s" 1 ). We report the 95% confidence upper- 
limits on emission lines from helium-like and hydrogenic Ne, 
Mg, and Si in Table 1 (upper limits on O emission lines are not 
constraining due to poor statistics at low energies). 

4. DISCUSSION 

We have analyzed the first CCD- and grating-resolution X- 
ray spectra of a millisecond X-ray pulsar in outburst. We find 
no convincing evidence for broad or narrow emission or absorp- 
tion features in the the EPIC-pn and RGS spectra (see Figure 1 
and Table 1). With our simple blackbody plus power-law model 
for the EPIC-pn spectrum, we find that the power-law compo- 
nent comprises 83% of the unabsorbed flux in the 0.5-10.0 keV 
band. The power-law component is hard: T= 1.44 ±0.01, 
but within a range typical for X-ray pulsars (e.g., T = 1.0—1.5; 
White, Swank, & Holt 1983). The EPIC-pn spectrum requires 
a blackbody with kT = 1 .05 ± 0.01 keV and an implied emitting 
radius of R = f x 3.01^ (rf/lOkpc) km. It is fair to say that 
our knowledge of the spectral hardening factor / is rather poor 
(for a discussion, see Lewin, van Paradijs, & Taam 1993). If / 
is less than 1.3, it would indicate that only part of the neutron 
star surface (a "hot spot") is involved. 

Markwardt et al. (2002) fit a simple absorbed power-law 
model to the RXTE spectrum of XTE J175 1-305, and mea- 
sure a softer power-law index (1.7 < T < 1.9). The effec- 
tive lower sensitivity bound of the RXTE Proportional Counter 
Array (PCA) is 3 keV; fitting the 3-10 keV EPIC-pn spec- 
trum with only a power-law, we find T = 1 .72 ± 0.01 for Nh = 
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(9.8 ±0.1) x 10 21 atoms cm -2 and T = 1.90 ±0.02 allowing N H 
to float (in this case, N H = [2.5 ±0.1] x 10 22 atoms cm" 2 ). Thus, 
the discrepant indices are easily explained in terms of different 
instrumental ranges. 

In modeling simultaneous Chandra/LETGS and RXTE spec- 
tra of XTE J0929-314, Juett, Galloway, & Chakrabarty (2002) 
measure similarly discrepant power-law indices. They suggest 
an astrophysical origin for the discrepancy in XTE J0929-314 
and XTE J1751— 305 in the form of a power-law with a break 
or a roll-over in the 1.4—4.4 keV range. A model consist- 
ing of blackbody and broken power-law components {Eh rea k = 
3.7 ±0.1 keVx, r £<3 . 7 = 1.34 ±0.03, r £>3 . 7 = 1.60 ±0.04) 
provides an improved fit to the spectrum of XTE J 175 1-305 
\x 2 /v= 1.090, v= 1882). Although these results suggest that an 
astrophysical origin for the apparent spectral evolution is pos- 
sible, a number of concerns remain. Given that: (1) the effec- 
tive area (flux) calibration of the EPIC-pn has 5% uncertainties 
(Kirsch 2002), (2) the LETGS (plus ACIS-S array) flux cali- 
bration has 10% uncertainties above 1 keV (H. Marshall, priv. 
comm.), (3) that RXTE PCA and High Energy X-Ray Timing 
Experiment (HEXTE) observations of the "Crab" nebula (taken 
to be a pure power-law for calibration) differ in power-law in- 
dex by S(T) ~ 0.1 (see, e.g., Wilms et al. 1999), and (4) that 
Markwardt et al. (2002) did not fit a blackbody component to 
the RXTE spectra of XTE J175 1—305, we do not consider an 
astrophysical origin for the discrepant indices to be required by 
the present data. 

The upper limits on the strength of narrow (consistent with 
instrument resolution; FWHM ~ 0.1 keV) and broad (FHWM = 
0.7 keV) Fe Ka emission lines are 4 eV and 6 eV, respectively 
(95% confidence). This stands in contrast the detection of a 
weak Fe Ka line in outburst spectra of the 401 Hz millisecond 
X-ray pulsar SAX J1808.4-3658 obtained with RXTE (Heindl 
& Smith 1998; Gierlinski, Done, & Barret 2002). It is possible 
that a strong, hot corona (or the magnetosphere, if the coranae 
in these LMXBs is small or disrupted) which produces to the 
power-law component may ionize the disk in XTE J 175 1-305 
to a degree which prevents Fe Ka emission. 

The absence of H-like and He-like resonance emission lines 
like those observed in 4U 1626-67 (Schulz et al. 2001) in the 



spectrum of XTE J1751— 305 is also consistent with a highly 
ionized disk in this source. Secondarily, the absence of such 
lines may be due to differences in the accretion flow geome- 
try due to the neutron star magnetic field. Whereas the mag- 
netic field in XTE J175 1—305 is likely to be similar to that in 
SAX J1808.6-3658 (B = (2-6) x 10 8 G; Wijnands & van der 
Klis 1998a), the magnetic field in 4U 1626-67 is likely much 
higher (B = 3x 10 12 G, Orlandini et al. 1998). The stronger 
magnetic field in 4U 1626-67 may disrupt the inner disk, allow- 
ing cool material in the outer disk to be irradiated by a weaker 
corona. Emission lines have also been seen in neutron star sys- 
tems which are viewed edge-on (so-called accretion disk coro- 
nae or "ADC" sources; see, e.g., Kallman et al. 2002); the 
absence of emission lines in XTE J175 1-305 may also be par- 
tially due to a low system inclination (Markwardt et al. 2002 
report no evidence of dips or eclipses). 

Studies of some ultra-compact systems suggest a Ne-rich 
companion (for observational results, see: Schulz et al. 2001; 
Juett, Psaltis, & Chakrabarty 2001; for theoretical discussions 
see Yungelson, Nelemans, & van den Heuvel 2002; Bildsten 
2002). Fits to the Ne photoelectric absorption edge in the RGS 
spectra of XTE J175 1-305 constrain the abundance of Ne to 
be at most 77% of the solar value along this line of sight (95% 
confidence upper-limit; see Section 3.3). This may suggest that 
a Ne-rich companion is unlikely in the case of XTE J 175 1-305. 
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Table 1 

Emission Line Upper-Limits 



Helium-like Hydrogenic 
Broad" Narrow* Broad" Narrow* 

(l(r 5 ph cm- 2 s~') (l(r 5 ph cm' 2 s" 1 ) (10~ 5 ph crrT 2 s"') (1(T 5 ph cnT 2 s" 1 ) 
Ne 1.6 1.4 0.9 0.5 

Mg 1.4 1.3 11 8.8 

Si 29 26 3.7 3.5 



NOTE. — The 95% confidence upper limits on highly ionized resonance emission lines within 0.1A of their expected rest-frame wavelengths in the RGS spectra 
of XTE J1751-305. " Assuming FWHM = 3000 km/s, similar to the Doppler-shifted pairs in 4U 1626-67 (Schulz et al. 2001). * The upper limits were calculated 
using Gaussians of zero width to model features narrower than the RGS resolution. Limits on the strength of O lines are not constraining due to poor statistics at low 
energies. These upper limits are similar to or below the fluxes measured from the lines detected in 4U 1626-67, indicating a constraining sensitivity was achieved in 
this XMM-Newton observation. 



